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ABSTRACT 

Waves and oscillations can provide vital information about the internal structure of waveguides 
they propagate in. Here, we analytically investigate the effects of density and magnetic stratification 
on linear longitudinal magnetohydrodynamic (MHD) waves. The focus of this paper is to study the 
eigenmodes of these oscillations. It is our specific aim to understand what happens to these MHD 
waves generated in flux tubes with non-constant (e.g., expanding or magnetic bottle) cross-sectional 
area and density variations. The governing equation of the longitudinal mode is derived and solved 
analytically and numerically. In particular, the limit of the thin flux tube approximation is examined. 
The general solution describing the slow longitudinal MHD waves in an expanding magnetic flux tube 
with constant density is found. Longitudinal MHD waves in density stratified loops with constant 
magnetic field are also analyzed. From analytical solutions, the frequency ratio of the first overtone and 
fundamental mode is investigated in stratified waveguides. For small expansion, a linear dependence 
between the frequency ratio and the expansion factor is found. From numerical calculations it was 
found that the frequency ratio strongly depends on the density profile chosen and, in general, the 
numerical results are in agreement with the analytical results. The relevance of these results for solar 
magneto-seismology is discussed. 

Subject headings: magnetohydrodynamics (MHD) - Sun: chromosphere - Sun: corona - Sun: oscilla- 
tions - waves 



1. INTRODUCTION 

High-resolution imaging and spectroscopic data from 
the Transitional Region and Coronal Explorer {TRACE) , 
Solar and Heliospheric Observatory (SOHO), Solar Ter- 
restrial Relations Observatories (STEREO) and Hinode 
have revealed a variety of wave modes in solar magnetic 
structures in the solar atmosphere (see, e.g., Banerjee et 
al. 2007; Aschwanden 2009; Taroyan & Erdelyi 2009; 
De Moortel 2009; Jess et al. 2009; Mathioudakis et 
al. 2011). Observations of standing slow magnetohy- 
drodynamic (MHD) waves have been reported by differ- 
ent authors, using SOHO/SUMER (Kliem et al. 2002; 
Wang et al. 2002, 2003, 2005), Yohkoh (Mariska 2005, 
2006) and recently, Hinode/ElS (Erdelyi & Taroyan 
2008). For an extensive review of the observation and 
modeling of standing slow magnetoacoustic waves see, 
e.g., Wang (2011). Coronal seismology, originally sug- 
gested by Uchida (1970), Zaitsev & Stepanov (1983), and 
Roberts et al. (1984), allows to obtain various physical 
parameters (e.g., magnetic field and density scale height) 
through matching the MHD theory and waves observa- 
tions in the corona. The concept was proposed to be used 
in any magnetic structure of the Sun by Erdelyi (2006) 
and labeled as solar magneto-seismology. The topic was 
extensively reviewed with plenty of references by Andries 
et al. (2009) and Ruderman & Erdelyi (2009). 

The theory of MHD wave propagation in magnetic 
structures in the solar corona has been developed by 
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modeling the magnetic structures as homogenous cylin- 
drical magnetic tubes embedded within a magnetic envi- 
ronment (Rae & Roberts 1982; Edwin & Roberts 1983; 
Roberts et al. 1984). Erdelyi & Fedun (2010) generalized 
these analytic efforts for fully compressible twisted mag- 
netic flux tubes. Perhaps the simplest model to study the 
oscillatory properties is a cylindrical flux tube in pressure 
balance, without the complexities of gravity, curvature, 
radiation, bulk motion, or heating. In this case, such 
modeling leads to a profile with constant pressure, den- 
sity, and temperature along the tube. 

However, more advanced equilibrium models have also 
been proposed with, e.g., dissipative effects and gravity 
(Mendoza-Briceno et al. 2004), and it has been found 
that the decay times of the standing slow modes are re- 
duced by 10% - 20% due to gravity when compared to the 
non-stratified loop models. Mendoza-Briceno & Luna- 
Cardozo (2006) included radiative cooling and heating on 
the study of longitudinal oscillations in hot, isothermal 
coronal loops with constant coronal heating; it was found 
that the lack of balance between cooling and heating does 
not affect the shape and decay time of the oscillations. 

Sigalotti et al. (2007, 2008) studied standing slow 
waves in hot coronal loops finding that in order to achieve 
the same rate of damping time as detected in the obser- 
vations, compressive viscosity has to be considered along 
with thermal conduction. Non-isothermal profiles of lon- 
gitudinal waves in hot coronal loops were examined nu- 
merically, arriving at longer periods and shorter damp- 
iyg^imes when the loop becomes more non-isothermal 
(^Erdelyi et al. 2008; Luna-Cardozo et al. 2008). 

However, the theory used in these previous studies as- 
sumed a constant magnetic field along loops. Verth & 
Erdelyi (2008) investigated the combined effects of mag- 
netic and density stratification on transversal coronal 
loop oscillations. It was found that even a relatively 
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small coronal loop expansion can have a significant and 
pronounced effect on the accuracy of the plasma density 
scale height measurements derived from observations of 
loop oscillations (see, e.g., Verth et al. 2008 for observa- 
tional case study). Soon after, Ruderman et al. (2008) 
applied the theory and found that the estimated coronal 
scale height is a monotonically decreasing function of the 
tube expansion, while studying transverse oscillations in 
a coronal loop with variable circular cross-sectional area 
and plasma density in the longitudinal direction. 

On average, the magnetic field strength is expected to 
decrease with height above the photosphere, although it 
has been very difficult to measure directly the variation 
of it in the corona. However, Lin et al. (2004) had made 
some progress using spectropolarimetry, and the results 
seem to confirm the decreasing with height of the mag- 
netic field strength. It is expected that most loops should 
expand with height above the photosphere, since the flux 
tube cross-sectional area and magnetic field strength are 
inversely proportional. This expansion is defined here by 



(1) 



where r a is the radius at the apex and rf is the radius 
at the footpoint (Klimchuk 2000). The loop expansion 
has been estimated for a number of loops, for example, 
Watko & Klimchuk (2000) reported mean values of 1.16 
and 1.20 for nonflare and postflare EUV loops, respec- 
tively, analyzing TRACE data. Klimchuk (2000) mea- 
sured a median value of T « 1.30 for soft X-ray loops 
using Yohkoh data. However, potential and magnetic 
field extrapolation had given larger loop expansions than 
the observed ones. DeForest (2007) suggested a possi- 
ble explanation for this based on the fact that resolu- 
tions of images have not been sufficient to actually detect 
coronal loop expansion. Regarding the chromosphere, it 
has been suggested theoretically that flux tubes must 
undergo significant expansion with height, the so-called 
magnetic canopy model, e.g., Gabriel (1976). However, 
to date there is little observational evidence to support 
this (Zhang & Zhang 2000). More recently, using the So- 
lar Optical Telescope (SOT) on board Hinode, for the 
first time Tsuneta et al. (2008) estimated an upper 
bound for chromospheric area expansion in the Sun's 
south polar region to be a factor of 345, giving a maxi- 
mum expansion factor of approximately 19 for chromo- 
spheric flux tubes. 

In this paper, the governing equation of the longitu- 
dinal MHD mode is derived and solved for two repre- 
sentative cases modeling solar atmospheric flux tubes: 
an expanding magnetic flux tube with arbitrary longitu- 
dinal plasma density and a density stratified flux tube 
with constant magnetic field. We examine the govern- 
ing wave equation within the limit of the thin flux tube 
approximation. The slow mode is decoupled from the 
other MHD modes, in a similar way as applied by Diaz 
& Roberts (2006), where the slow mode was studied in 
density-structured coronal loops with constant magnetic 
field. The purpose of this paper is to quantify the sepa- 
rate effects of the expansion factor T and density stratifi- 
cation. The shooting method is applied to find numerical 
solutions of the general wave governing equation to com- 
pare them with analytical approximations. 



2. MAGNETIC FIELD EQUILIBRIUM CONFIGURATION 

The magnetic field equilibrium that decreases in 
strength with height above the photosphere is modeled 
by an expanding flux tube with rotational symmetry 
about the z-axis in cylindrical coordinates (r, 9, z). Ne- 
glecting curvature along the tube axis we model an ex- 
panding tube with a straight central axis, i.e., a magnetic 
bottle. The tube ends are frozen in a dense photosphcric 
plasma at z = ±L, and the flux tube has an arbitrary 
density depending on z. This expanding magnetic field 
in equilibrium has two components 



B = B r {r, z)e r + B z (r, z)e z 



(2) 



Following the derivation of Verth & Erdelyi (2008) for 
a potential field configuration, if small expansion is as- 
sumed it is possible to obtain an explicit expression for 
the perpendicular distance from the tube axis to a mag- 
netic surface with footpoint distance from the axis de- 
fined by rf, i.e., 



r (z) « r f <^ 1 + 



(1 - T 2 ) [cosh(z/L) - cosh(l)] 



r 2 



1 — cosh(l) 



-1/2 



(3) 



Note that in the thin tube approximation Equation ([3]) is 
only a function of z. Similarly, the magnetic field com- 
ponents B r and B z at near the tube axis can also be 
described explicitly as a function of z, 



B r {z) « -B z .f 



(1 -T 2 ) sinh(z/L) r (z) 



2r 2 



1 - cosh(l) L 



(4) 



and 

B z (z) 



B 7 



(1 -T 2 ) [cosh(z/L) 



r 2 



3Sh(l)] ] 

1 - cosh(l) J 



Therefore, B r and B z are related by, 



1 , .dB z 
B r (z)*~-r (z) — 



(5) 
(6) 



3. GOVERNING EQUATIONS 



The ideal MHD equations are linearized by assuming 
small magnetic perturbations b = (b r ,Q,b z ) and velocity 
perturbations v = (v r ,0,v z ) about a plasma in static 
equilibrium. In the derivation we neglect gravity and as- 
sume constant kinetic plasma pressure. Note this means 
along our model of solar flux tubes that the plasma is 
not isothermal and the assumption of constant plasma 
pressure has greater validity in the corona than in the 
chromosphere. 

Following, e.g., Roberts & Webb (1978), Roberts 
(2006), and Diaz & Roberts (2006), we also neglect the 
effect of the external environment on the perturbations, 
i.e., we assume the tube is in a quiescent environment. 
This means we do not consider any external forces act- 
ing on the tube. Inclusion of such effects is essential in 
the studies of, e.g., p-mode absorption of photospheric 
flux tubes (see Bogdan et al. 1996). The resulting MHD 
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(7) 

(8) 

(9) 
(10) 

(11) 
(12) 

(13) 



where t is time, r and z are the radial and longitudinal 
coordinates in the tube, po is the plasma mass density in 
equilibrium, Pt is the total perturbation to pressure, p is 
the kinetic pressure perturbation, po is the kinetic plasma 
pressure in equilibrium, B r and B z are the background 
components of the magnetic field, 7 is the ratio of specific 
heats, and p is the magnetic permeability. 

3.1. Magnetic Flux Tube Equilibrium 

The potential magnetic field configuration chosen in 
Section [2] and our choice of constant plasma pressure put 
restrictions on the possible types of flux tube equilibria 
we can model. It is always assumed that the external 
magnetic field is balanced by the internal one, therefore 
the flux tube models are not in a magnetic field free en- 
vironment. Then, our assumption of constant plasma 
pressure demands that a flux tube with greater internal 
plasma density than external one must also be cooler 
than its environment. On the other hand, if the tube is 
less dense than its environment, the internal temperature 
must be hotter than the external. 

3.2. Velocity Wave Equation 

Implementing the thin tube approximation we intro- 
duce the spatial and temporal scalings r = ei?, z = Z and 
t = t, where e< 1. We show the details of our deriva- 
tion for the governing thin tube v z velocity equation in 
Appendix El Let us introduce the following Fourier de- 
composition v z (Z,t) = v z {Z) exp(— iwr), where ui is the 
angular frequency of the oscillations. Then, the second- 
order ordinary differential equation for the longitudinal 
velocity wave is 



alZ 2 



where 



h{z) 



■3- + MZ) 



= 0, 



(14) 



1 dB z 
B z dz 



h{Z) 



1 d 2 B z 
B z dZ 2 



1 (dB z 



m V dz 



For constant magnetic field, Equation (fT4)) reduces to 
d 2 v z [uj 2 \ 

^+Uh =0 ' (15) 

for 



— = w po 



_p_ 



1 



(16) 



where the density may depend on Z. This equation (| 15[) 
agrees with the equation obtained by Diaz & Roberts 
(2006), where B z was considered constant. The solutions 
to Equation (|15p will depend on the functional form cho- 
sen for po in Equation (TT6|) . 

3.2.1. Slow Modes in a Homogeneous Tube 

For checking the derivations, we may recover the find- 
ings for a homogeneous tube (e.g., Edwin & Roberts 
1983). In a straight magnetic flux tube with con- 
stant density, the solutions to Equation (|15[) are simply 
trigonometric functions. Applying the line-tying bound- 
ary at the ends of the tube, v z (±L) = 0, we find the 
frequencies of the even modes given by 



ct (2b — 1)tt 



n = 1,2,.. 



L 2 

whereas the odd modes have frequencies given by 

o C T 



L 



nn, n = 1,2, 



(17) 



(18) 



The ratio of frequencies of the first overtone and funda- 
mental mode is equal to 2, as expected. 

4. ANALYTICAL SOLUTION FOR A SMOOTH DENSITY 
PROFILE 

Since we are using cylindrical coordinates we model a 
solar coronal loop by a straight axisymmctric magnetic 
flux tube. Therefore we neglect the effect of flux tube 
curvature and are simply modeling the effect of gravita- 
tionally stratified plasma in a coronal loop which would 
produce a density profile symmetric about the loop apex. 
The tube length is 2L and its radius is r . The plasma is 
permeated by a uniform magnetic field B directed along 
the tube axis, B = B z z. The plasma density, po(Z), is 
greater at the loop footpoints than at the apex so this is 
approximated by the function 



Po(Z) 



Pa 



[l-(l- Kl )(Z/L) 2 } 



2 • 



(19) 



where p a = po(0) and pt — po(±L) are the apex and 
footpoint densities, respectively, and «i = (pa/pt) 1 ^ 2 the 
stratification parameter. The solid line in Figure Q] shows 
this smooth density profile as a function of Z . In Section 
we choose a more applicable and realistic exponen- 
tial density profile and solve Equation (1151) numerically. 
However, the choice of density given by Equation (TT9"]) al- 
lows us to have a straightforward analytical insight into 
the effect of density stratification on longitudinal oscilla- 
tions. 



4 



Luna-Cardozo, Verth, & Erdelyi 




o 

-L L 

z 

Fig. 1. — Density pr ofile of the loop. The solid line shows the 
smooth density profile H91 while the das hed line shows the expo- 
nential density profile of a coronal loop i]44tl . and the dotted line 
show s the exponential density profile of a chromospheric flux tube 
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Fig. 2. — Dependence of the frequency ratio of the first over- 
tone and fund am ental mode on Ki, for a semi-circular coronal loop 
in solid line 1^251 1, and for a vertical chromospheric flux tube in 
dashed line 1 140 f . The "+" and "x" symbols represent the numer- 
ical solutions of Equation H20H for Z £ [—L,L] and Z £ [— L, 0], 
respectively. 



With the equilibrium density profile given by Equation 
T51 ). Equation (|T5|) takes the form 



1 



dZ 2 4, a [l-(l-/ei)(Z/£) 2 r 



0, 



(20) 



for cx, a = ct(0), the tube speed at the apex. The general 
solution of this equation is (see Polyanin & Zaitsev 2003) 



v z = v/l-(l- Kl )(Z/L) 2 (Ci cos(w) + C 2 sin(u)) . 



where C\ and Ci are arbitrary constants, and 



(21) 




To study a standing wave the boundary conditions 
v z (±L) = are applied, and are satisfied when either 
Ci = and cos[u(L)] = 0, or C\ — and sin[u(L)] = 0. 
The first condition corresponds to even modes, and the 
second to odd modes. The frequencies of the even modes 
are given by 



(2n- lJV In 



K) 2 = 
1 + VT 



(f - K X ) C \ . 

L 2 



f-x/T - ^ 



(23) 



for n — 1,2,..., while the odd modes have frequencies 
given by 



(2mrY In 



w? = — 

1 + VI - hi 



Kl)c| . 



L 2 

-2 



1 - VT^T 



(24) 



for n = 1,2, ... 

The frequencies of the fundamental mode and the first 
overtone are given by Equations (|23p and (IM1) with 
n = 1, respectively. Theoretically it is predicted that 
the frequencies of higher harmonics have a much stronger 
dependence on density stratification (e.g., Andries et al. 



2005). The ratio of frequencies of the first overtone and 
fundamental mode is given by 



L0 2 



4tt 2 



In i+^I 
1— VI— «i 



IT 2 + 



1/2 



(25) 



which is the same result obtained for the frequency ratio 
of the transversal mode by Dymova & Ruderman (2006). 
This is a remarkable property of MHD oscillations in 
structured waveguides. This is also rather assuring, as 
we arrived to these results using a completely different 
approach and modeling. The dependence of this ratio of 
frequencies on n\ is shown in Figure [5] by the solid line 
(semi-circular coronal loop). If we consider the limit of 
a non-stratified loop, i.e., k\ — ¥ 1, we find the frequency 
ratio tending to 2, and it can be approximated by 



^ = 2-4(1-*). 



(26) 



Equation ([26|l shows that the frequency ratio ui-xjuix < 2 
for a stratified loop with constant magnetic field. 

5. FLUX TUBE EXPANSION WITH CONSTANT DENSITY 

Let us now study the effect of magnetic stratification 
with a constant density on the longitudinal oscillations. 
For this configuration it is more convenient to use the 
governing pressure perturbation equation. See Appendix 
[B] for insight into the relationship between v z (Z) and 
p(Z) when there is l ong itudinal stratification. The pres- 
sure wave equation (|B1|) with constant density is 



d 2 p 
dZ 2 



1 dB z dp 
~K~dZdZ 



0. 



(27) 



We can transform Equation (|2T|) to its canonical form 
using the change of variable p = P\J B z . Then, Equation 
(l27l) becomes 



d 2 P 
dZ 2 ' 



uj 2 1 1 d 2 B z 

~72 



2 B r dZ 2 



3 1 

1b 2 



dB z 
dZ 



P = 0. (28) 



Oscillations in solar waveguides 



5 



For weak magnetic stratification the uj 2 jc\ term has a 
dominant effect on the eigenvalues. To analytically in- 
vestigate the behavior of the eigenvalues in this regime, 
we approximate Equation (|28|) with 



(29) 



dz 2 + 4 



We suggest to use a rational function for the tube speed 
ct defined by 



Ct(Z) 



C T,i 



CT,{ 
CT.a 



1 2 



(30) 



for CT, a = ct(0) and cr,f = ct(±£) being the apex and 
footpoint tube speeds, respectively. In Section[71 we solve 
the governing velocity equation (fT4")l with our potential 
field definition of B z (z) from Equation ([5]) numerically. 
However, the choice of tube speed given by Equation 
(|3"D|) allows us to have a straightforward analytical insight 
into the effect of magnetic stratification on longitudinal 
oscillations. 

An exact solution to Equation (|29p with ct defined by 
Equation (|30|) is given by (see Polyanin & Zaitsev 2003) 



P = yJa 2 Z 2 + I (C 3 cos(zy) + C 4 sin(z/)) , 
where C3 and C4 are arbitrary constants, with 



2„2 



arc 



I arctan(aZ), 



and 



CT,f 
CT,a 



1. 



(31) 



(32) 



(33) 



In this case, to find the standing mode solution the same 
boundary conditions v z (±L) = P'(±L) = are applied. 
The frequencies of the even modes are given by 



= n; 

CT,a 

for n = 1,2, .. 
are 



or,. 



+0(aL) 4 



T + i) (ai)2 + 0(ai)4 ' (34) 

while the frequencies of the odd modes 
(n — |)tt 1 



2(n-|)7T 



(<xL) 2 
(35) 



for n = 1, 2, .... The ratio of frequencies of the first over- 
tone and fundamental mode is 



UJ-2 



r (aL) 2 +0(aL) 4 



(36) 



We now want to see how the effect of T is related to 
plasma ft. Since we have not defined T explicitly in Equa- 
tion pop , we combine this definition of tube speed with 
our potential magnetic field of B z {z) given by Equation 
([5]), which results in the ratio of tube speeds, 



CT,f 
CT,a 



2 + 7&r 4 

2 + 7/3f 



1/2 



(37) 
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(a) Semi-circular coronal loop 
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(b) Vertical chromospheric flux tube 
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Fig. 3. — Dependence of the ratio of freq uenc ies of t he first over- 
tone and fundamental mode (Equations J38D and 1431 ) on the 
expansion parameter T for different values of /3{. Solid, dashed, 
dotted, dot-dashed, triple-dot-dashed, and long-dashed lines corre- 
spond to /3f = 0.1, 0.5, 1, 2, 5, and 10, respectively. 

where fti = 2[ipo/B 2 { is the beta plasma at the foot- 
points. Therefore, Equation (131)]) shows that change in 
frequency ratio is dependent only on the loop expansion 
factor r for any finite beta plasma. 

If we consider the limit of a non-expanding loop, i.e., 
r lj we recover a constant tube speed along the loop. 
The ratio of frequencies can be Taylor expanded for r ~ 1 
and 7 = 5/3, giving 



^ = 2 



30ftt 



(6 + 5/3 f )7T 



-(r-i) 



(38) 



Equation (|38p clearly shows that uj 2 /uj\ < 2 for oscilla- 
tions in an expanding magnetic flux tube with constant 
density, and it shows a clearly linear dependence with the 
expansion factor. This relationship between ujz/uix and 
r is shown in Figure |3U a). It is clear that increasing the 
magnetic stratification takes to a lower frequency ratio, 
and this effect is more significant for plasmas with higher 
fti- 

6. APPLICATION TO SOLAR PHYSICS 

Our results are relevant to magneto-seismology, e.g., 
estimating the coronal density scale height by using the 
observed ratio of the fundamental frequency and first 
overtone of longitudinal loop oscillations. 

6.1. Application to the Corona 
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In the solar corona the thermal pressure is generally 
smaller than the magnetic pressure, giving a plasma /? 
parameter <C 1. Stratified coronal loops (0 < «i < 1) 
with constant magnetic field give a frequency ratio lower 
than 2 (Equation flU). 

For a 1 MK average corona, the expected hydrostatic 
scale height H should be about 50 Mm, therefore for 
any coronal loop with length of the order of H or less 
could be approximated by the configuration in Section 
5, where the density is constant inside the loop and the 
tube speed (i.e., geometry of the magnetic field) is given 
by Equation ([50]). 

For a very small beta plasma (i.e., 0.01 < j3 < 0.1) 
the expansion has a weak effect on the frequency ratio, 
as it can be seen in Figure [3ja) for the solid line. The 
frequency ratio can be Taylor expanded for a very small 
plasma /3f parameter (/3f <C 1), giving 



LU-2 
UJl 



5 fT A 

4 



1 



7T 



(39) 



Equation (I39[) shows that LO2/UJ1 < 2 for an expanding 
magnetic flux tube (r > 1) in the corona with constant 
density. 

6.2. Lower Solar Atmosphere 

We now calculate the frequencies of standing modes 
in an expanding chromospheric (/3 > 1) flux tube with 
fixed boundaries at the photosphere and transition re- 
gion, since, e.g., Fujimura & Tsuneta (2009) have re- 
cently obtained observational evidence of such waves. 
Note that it could be possible that flux tubes undergo 
large expansions in the chromosphere (e.g., Tsuneta et 
al. 2008), so our assumption of weak magnetic strati- 
fication may have limitations in application to chromo- 
spheric wave observations. However, the size of correc- 
tions to eigenfrequencies if larger flux tube expansions 
are considered are as of yet still unquantified and must be 
the focus of a future study. Regarding our current model, 
to calculate the eigenfrequencies of a vertical chromo- 
spheric weakly expanding flux tube, we need to solve the 
eigenvalue problem in only half of our magnetic bottle, 
and therefore, the boundary conditions v z (—L) = and 
v z (0) = are applied. 

For the loop with density stratification and constant 
magnetic field, it is found that C\ = and sin[u(— L)\ = 
for any arbitrary constant Ci- The frequencies for all 
odd and even modes (w„) are given by Equation (|24[) . 
The frequencies of the fundamental mode and the first 
overtone are given by Equation (|24|) with n = 1 and 2, 
respectively. The ratio of frequencies of the first overtone 
and fundamental mode is 



ijJ-2 



-in- 



1/2 



(40) 



It is important to note, that the value of uji/uji strongly 
depends on which functional form is chosen for the equi- 
librium density. Figure [5] shows a comparison of the 
dependence of the frequency ratios on the parameter 
Ki, in the lower atmosphere (dashed line) and in the 
corona (solid line) . When we consider the limit of a non- 
stratified loop, i.e., Ki = 1, we recover the ratio of fre- 



quencies equal to 2. Equation (1401) can be approximated 
to 

- = 2-^(1-^ (41) 



for Ki ~ 1. Equation (|41l) shows that the frequency 
ratio UJ2/0J1 < 2 for a stratified loop (k\ < 1) in the 
lower solar atmosphere. Note that for our coronal loop 
model the densities at both tube ends are equal with 
the density profile symmetric about the apex. In our 
vertical chromospheric flux tube model the densities at 
both ends are not equal, and the density is monotonically 
decreasing as a function of height. It is found that the 
effect of the density stratification on the frequency ratio 
is larger in the corona than in the lower atmosphere, due 
to the asymmetric nature of the density profile in the 
chromosphere. 

In the second case, a loop with magnetic stratification 
and constant density, the same boundary conditions are 
applied. The frequencies for all odd and even modes (u> n ) 
are given by Equation (|34|) . 

The ratio of frequencies of the first overtone and fun- 
damental mode is 



^ = 2--L(aL) 2 + 0(aL) 4 . 
The Taylor expansion of Equation (|42)) for V ~ 1 is 



UJi 



15 



ft 



2 (6 + 5/3 f )7T 



-(r-i), 



(42) 



(43) 



showing again that UJ2/UJ1 < 2 and it has a linear de- 
pendence on r for standing waves in the lower atmo- 
sphere with constant density and weak magnetic expan- 
sion. This frequency ratio as a function of T is shown 
in Figure EJb) . The effect of the magnetic stratification 
on the frequency ratio is smaller in the lower solar atmo- 
sphere than in the corona. 

7. NUMERICAL SOLUTIONS 

In this section, we compare the analytical approximate 
solutions with the numerical solution of Equation (I14p . 
using the shooting method based on the Runge-Kutta 
technique. 

7.1. Stratified Loop with Constant Magnetic Field 

Equation (fbfj) becomes Equation (|15l) when B z is con- 
stant. Equation (fT5")) can be solved numerically for both 
a semi-circular coronal loop and a vertical chromospheric 
flux tube, depending on the equilibrium chosen. 

In the first case, for a semi-circular coronal loop per- 
pendicular to the plane of the photosphere with z 6 
[-L, L], the density is defined by (see dashed line of Fig- 
ure [T]) 



Po(z) = Pfexp 



2L Ckz\ 

cos — 

ttH \2LJ 



(44) 



where H is the density scale height and pf the footpoint 
density. Hence, in this case the parameter ki of Section 
4 is a function of H and is given by 



exp 



2L 
"irH 



(45) 



The frequency ratio obtained with the density profile (|44l) 
as a function of Ki is shown in Figure [4] by the solid line. 
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Fig. 4. — Dependence of the ratio of frequencies of the first over- 
tone and fund ame ntal mode on ki, for a semi-circular coronal loop 
in solid line q44t and for a vertical chromospheric flux tube in 
dashed line <46t . obtained from the numerical calculations. 

It is clearly lower than 2, similar to the results obtained 
in the analytical case. 

For a vertical chromospheric flux tube with z G [-L, 0], 
the density profile shown by the dotted line in Figure [1] 
is 

(z + LY 
Pfexp — — 



Po(z) 



(46) 



where the parameter K\ is now given by 



L 



exp [ -- 



(47) 



Equation (|15l) with density profile given by Equation (|46l) 
has the well-known solution 



v z {Z) = C 5 J 



/ 2Hlu 



+ C 6 Y 



2Hu 
ct{Z) 



(48) 



where C5 and Cq are arbitrary constants. To find the 
eigenvalues using solution (|48|) , we must solve a transcen- 
dental equation by either analytical or numerical tech- 
niques. Here, we choose to solve it numerically. See, e.g., 
McEwan et al. (2008) and Verth et al. (2010) for analyt- 
ical solutions of equivalent equations in their studies of 
other MHD waves in limiting cases of weak stratification 
for both coronal loop and vertical flux tube geometries. 

The dependence of the ratio of frequencies on k,\ for the 
density profile (|46|) is shown in Figure [4] by the dashed 
line. In this case, the frequency ratio is greater than 2, 
confirming that this parameter strongly depends on the 
choice of the functional form of density. This means that 
caution must be used when interpreting the frequency ra- 
tio of chromospheric standing modes. For example, the 
choice of a density profile that gives a tube speed increas- 
ing linearly with height results in 1^2/ uji < 2, while the 
choice of a density profile giving a tube speed exponen- 
tially increasing with height gives CJ2/W1 > 2. 

7.2. Expanding Loop with Constant Density 
Using Equation ([5]) for B z (z) with z € [— L, L] we can 



now compute the numerical solution of Equation (|14|) for 
longitudinal oscillations in a coronal loop. Figure [Ha) 
shows the frequency ratio of the first overtone and fun- 
damental mode as a function of the expansion parameter 
r for different values of /3f. It is found that increasing the 
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(b) Vertical chromospheric flux tube 
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Fig. 5. — Dependence of the frequency ratio of the first overtone 
and fundamental mode on the expansion parameter T for differ- 
ent values of /3f obtained from the numerical calculations. Solid, 
dashed, dotted, dot-dashed, triple-dot-dashed, and long-dashed 
lines correspond to /3f = 0.1, 0.5, 1, 2, 5, and 10, respectively. 

magnetic stratification leads to a lower frequency ratio, 
and this effect is more significant for solar waveguides 
with higher /3f . These results are very similar to the re- 
sults obtained in the analytical analysis. 

To study the eigenmodes of a vertical chromospheric 
flux tube, we use Equation ([5]) for B z (z) with z e [— L, 0] 
obtaining the numerical solution of Equation (|14l) . The 
dependence of the frequency ratio of the first overtone 
and fundamental mode on the expansion parameter is 
shown in Figure [5jb). The effects of the expansion T 
is somewhat smaller in the chromosphere than in the 
corona, being consistent with the analytical results. Also, 
the values of the frequency ratio for T < 1.5 (weak expan- 
sion) obtained from numeric and analytical calculations 
are very similar. 

8. IMPLICATIONS FOR SOLAR MAGNETO-SEISMOLOGY 

The governing equation flT4"]) is valid for both density 
and magnetic stratification. We investigated the two ef- 
fects, density and magnetic variations, on the value of 
UJ2/0J1 separately for clarity in Sections!?] - 7. Now we 
discuss the combined effects in realistic solar waveguides 
where it is likely that the two types of stratification are 
present simultaneously. 

8.1. Semi- circular Coronal Loops: Effect of both 
Density and Magnetic Stratification 

In the case of a coronal loop, it was shown in Sections 
|4] and [5] analytically and Section [7] numerically, that if 
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there is density stratification (magnetic field constant) 
then W2/W1 < 2 and if there is magnetic stratification 
(constant density) then also UJ2/0J1 < 2. Hence, regard- 
ing the value of 0J2 /<-0i the two effects are not competing 
against each other, in contrast to the kink mode (see, 
e.g., Verth & Erdelyi 2008; Ruderman et al. 2008). This 
result is robust, in particular, in the functional forms 
chosen for po{Z) and B Z {Z) as long as there is symme- 
try about the loop apex. 

Observational claims of estimating slow mode values of 
from intensity perturbations in coronal loops have 
been made using both EUV (Srivastava & Dwivedi 2010) 
and X-ray (Kumar et al. 2011) data from Hinode. Sri- 
vastava & Dwivedi (2010) claim to detect two separate 
cases in EUV coronal loops of 0-12/^1 < 2. However, in 
these observed coronal loops the value of plasma f3 may 
be very low, and in this limit ct(z) ~ c s (Z). In this sit- 
uation it is the variation of c s {Z) in wave equation (|14|) 
that would have the dominating effect on the value of 
UJ2/1JJ1 and magnetic stratification, even if present, would 
not play a significant part. Assuming the observed loops 
are isobaric and density is decreasing as a function of 
height, then even in the zero-/? limit, our model gives 
OJ2/0J1 < 2, consistent with these particular EUV obser- 
vations. Physically, this would be due to c s (Z) increasing 
with height in a coronal loop, i.e., temperature hotter at 
the apex than at the footpoints. 

Regarding the study by Kumar et al. (2011) of slow 
mode values of uj2/lo\ in X-ray small-scale loops, Ta- 
ble 1 of that paper shows six separate estimates of 
and one particular case of uj^/lj\. Three esti- 
mates have LJ2/UJ1 < 2, similar to the measurements of 
Srivastava & Dwivedi (2010), but in contrast three also 
have > 2. In the low-/3 limit, the observations of 

CO2/0J1 > 2, could be explained by c s (Z) decreasing as 
a function of height in these loops, i.e., the temperature 
being hotter at foot-points than at the apex. However, 
to explain this our isobaric model would require den- 
sity increasing with height to have c s {Z) decreasing with 
height and this is unphysical. To model such a coro- 
nal loop more realistically, one would also need to model 
p(Z) and choose p(Z) such that temperature is hotter at 
foot-points than at apex. This would be an important 
future extension to this current work. 

8.2. Vertical Flux Tubes: Effect of both Density and 
Magnetic Stratification 

Our study of magnetic stratification (constant density) 
in vertical flux tubes resulted in ^2/^1 < 2 for both the 
hyperbolic profile for B Z (Z) given by Equation §5$ and 
the tube speed given by the rational function in Equa- 
tion (|30|) . However, it was shown that for the case of 
density stratification (constant magnetic field) in a ver- 
tical flux tube, LJ2/UJ1 could either be greater or less 
than 2 depending on the functional form chosen. E.g., 
W2/W1 < 2 for the smooth rational function of density 
given by Equation (| 19[) and ^2/^1 > 2 for the expo- 
nential profile of po(Z) given by Equation (|46|). In our 
investigation we have shown that, depending on the func- 
tional forms chosen, density stratification and magnetic 
stratification could both cause the value of W2/W1 to be 
less than 2, or the effects could oppose each other. This 
means that magneto-seismology for vertical flux tubes in 
the chromosphere may be more subtle than in the case of 



coronal loops when it comes to interpreting the observed 
values of 1^2/^1- 

9. SUMMARY AND CONCLUSIONS 

In this paper the effects of density and magnetic strat- 
ification on longitudinal oscillations in isobaric coronal 
and chromospheric conditions have been studied. So- 
lar waveguides were modeled as axisymmetric cylindrical 
magnetic tubes. The governing equation were derived 
and solved by analytical approximation, and examined 
in the thin flux tube limit. We studied the effects of 
magnetic stratification while the density is constant and 
density stratification for a constant magnetic field. 

From the analytic solutions, both density stratified and 
expanding coronal loops have 102/1^1 < 2. For small ex- 
pansion, a linear dependence between the frequency ra- 
tio and the expansion factor is found. It was also found 
that the effect of magnetic field strength decreasing with 
height has the same effect on the frequency ratios to that 
of gravitational density stratification, in contrast with 
the results for kink modes (Verth & Erdelyi 2008; An- 
dries et al. 2005). 

It was found that the introduction of waveguide struc- 
turing results in a modification to the oscillatory fre- 
quency of the mode. The expression for the frequency 
ratio obtained in Section 5 (Equation (|3E])) depends on 
the expansion parameter T and an additional dependence 
on f3i was found, due to the inclusion of kinetic pressure 
in our model. 

Next, numerical solutions were performed. The nu- 
merical results were consistent with the analytic solu- 
tions for the coronal loop model. Also, for the vertical 
chromospheric flux tube model it was found that the fre- 
quency ratio strongly depends on the functional form of 
the density, suggesting that caution must be used when 
interpreting the frequency ratios of chromospheric stand- 
ing modes. 

The effect of gravity was neglected in the present study. 
It is known for propagating longitudinal waves that the 
cutoff frequency is increased by the inclusion of gravity 
(see, e.g., Roberts & Webb 1978), which introduces a 
Brunt- Vaisala (buoyancy) term into the governing equa- 
tions. This term has been shown to be relevant for the 
leakage of p-mode driven longitudinal waves into the up- 
per atmosphere (De Pontieu et al. 2004, 2005), since 
the influence of the Brunt- Vaisala term in the cutoff fre- 
quency can be reduced by the amount of magnetic field 
inclination, i.e., the more tilted the field from the verti- 
cal, the lower the effective cut-off frequency. In a future 
study, it would therefore be of great interest to inves- 
tigate the effect of gravity and field inclination on the 
longitudinal standing modes, e.g., to quantify the im- 
portance of the Brunt- Vaisala (buoyancy) term on the 
eigenfunctions and eigenvalues. 

It is important to progress the field of magneto- 
seismology that any proposed model can be tested 
against the observed oscillatory properties of solar waveg- 
uides. The theory presented in this paper should be help- 
ful in this regard, modeling both magnetic and density 
stratification along solar waveguides which can be used to 
interpret observations of standing slow oscillations such 
as those by Srivastava & Dwivedi (2010) and Kumar et 
al. (2011). In the future, it will also be beneficial to find 
the solutions of the governing equation when stronger 
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magnetic stratification is considered. 
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APPENDIX 

DERIVATION OF VELOCITY WAVE EQUATION 

In this Appendix we derive the velocity wave equation (| 14|) . Equations (|7|)- (|13p may be combined to yield the general 
equations for v z , v r , and Pt, 



d_ (dPr\ 
~d~z \~aT) 



d 2 v z 



l>i> 



B 2 r d 2 v z B r B z d 2 v 



*2B r 0B Z 

fi dr 



dt 2 
B r d 2 B z 
ix drdz 

B r B z 
fir 



fi dr 2 fi 
+ 1 dB z dB r 
fi dr dz 



drdz 
B z d 2 B 



8v r 

dr 



' B r dB; 

fi dr 

*B 2 d 2 V 



2B Z dB 
fi dz 



B z dB, 
fi dr 
B r B z d 2 v 



2B Z dB : 
fi dz 
d fdP T \ 
Ik- V~9T/ 

dv r 

~dz~ 



B r dBj 
fi dr 

1 f dB 2 

fi dz 2 fi \ dz 

B r dB z B r d 2 B : 



B z dB r 
fi dz 



dv z 
dr 

B r B z d 2 v 



B z dB z ^ B r dB : 
fi dz fi dr 
B 2 d 2 v. 



fir dr 

d 2 V r 



po- 



+ B z d 2 B z 

fi dr 2 fi drdz 
B 2 Z d 2 v r B r B z d 2 v r 



dv z 
~dT 

B 2 Z dv r 

fir dz 



dt 2 fi dz 2 
1 dB z dB r B r d 2 B 
fi dr dz 



fi drdz 



2B Z dB. 
fi dz 
Bl 



dv z 
~dz~ 



2B r dB., 
fi dr 

dP T _ _B| f dvr 
dt fi \ dr 

B r dB r B z dB, 
fi dz fi dz 



dr 2 fi drdz 

B r B z 2B Z dB; 
fir 2 fir dz 

+ \B Z dB r _ B r dB; 

fi drdz] \_ fi dz fi dz 
B z d 2 B z B z dB r 1 dB r dB-. 

fir dz fi dr dz 
B z d 2 B T B r dBi 
fi dz 2 fir dz 

B r B z dv z B 2 dv 2 
fi dr fi dz 
B r B z dv r (B r dB z B z dB : 
fi dz \ fi dz fi dr 



= 0, 

dv r 

dr 



fi drdz fi dz 2 



B r d 2 B r 



fi drdz 
— IPO (V • v) + 



(Al) 



(A2) 



(A3) 



If we compare Equations (| Al|) - (| A3|) with Equations (l)-(3) from Diaz & Roberts (2006, who studied density stratifi- 
cation with a constant magnetic field) , the added complexity to the governing wave equations if we have an equilibrium 
with an expanding magnetic field can readily be appreciated. 

Here we are interested in the behavior of slow magnetoacoustic modes, we introduce the following spatial and 
temporal scalings: 



tR, z = Z, v r 



, B r 



eBf{, B z = B Zl and t = r. 



(A4) 



It is interesting to note that such scaling have been used in the analysis of slow modes within resonant layers in a 
magnetically structured plasma (e.g., Ballai et al. 1998) and to study slow MHD waves in a stratified medium (e.g., 
Roberts 2006, except that in his study the stretching was only on the position and time since the magnetic field was 
considered constant). 

Considering the thin tube approximation, e <C 1, i.e., the r-coordinate has a small range in comparison with the 
z-coordinate. Since we are interested in the longitudinal component of the wave, we focus on the equation of v z , and 
after the scaling it becomes 



Po 



d 2 v z 



c 2 s Bl d*v z 



+ 



4 2B R dB z 



dr 2 1 8Z 2 
B R dB R , B z dB R 



c\ [i OR 2 



H dZ c| \ fi dR fi dZ 



dv z 
~dR 



c 2 2B R B z d 2 v z 
cf fi dRdZ 

c 2 B z d 2 B z 







B z dB z - 


dv z 


[U f 2 


4 j 


li dZ _ 


dZ 



c 2 fi dZ 2 



dB z 
dZ 



v z = 0, (A5) 



i 1 /2 j 

where ca = (B 2 /pp ) and c s = (7Po//°o) 1/z are the Alfven and sound speeds, respectively. The square of the fast 
phase speed is defined by ef = c 2 + c\ and the tube speed is given by c^ 2 = c^ 2 + c^ 2 . 

The R- derivatives of perturbed quantities are very small compared with the Z-derivatives since we consider the thin 
tube approximation, i.e., dv z /dR <C 1. 



,1/2 



PRESSURE WAVE EQUATION 

In this Appendix, we present the governing pressure equation and the relationship between velocity and pressure 
amplitude for longitudinal stratification. Equations (|7))- (fT3|) may also be combined to give an equation for the (per- 
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turbation) pressure amplitude p(z), defined by p{z, t) — p{z) exp(iuit). If one follows this route, one will arrive at 



d 2 p 
dZ 2 



1 dpo 
Po dZ 



1 dB z 
B z dZ 



dp 
~dZ 



-p = 0, 



(Bl) 



after applying the same scaling as defined by (|A4|) . Equation (|B1|) is consistent with the results obtained by Roberts 
& Webb (1978). From Equation ((SJ), the following relation is obtained: 



dp , 7 , dv z 

We can find the velocity amplitude equation from the pressure equation by using (|B2 



(B2) 



d 2 v z 
dZ 2 



1 dc 2 T 



1 dp a 
Po dZ 



1 dB z 
B z dZ 



dv z 
~dZ 



1 d 2 B z 



B, dZ 2 



1 

m 



ob 2 



dZ 



1 d Po 1 dB z 



po dZ B z dZ 

which is indeed equivalent to Equation (| X4[) and that can be checked rather easily. 



1 dc\ 1 dB z 
^~dZB^~dZ~ 



v z =0, 



(B3) 
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